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ABSTRACT 

We consider the decoding of bit interleaved coded modulation (BICM) 
applied to multiband OFDM for practical scenarios where only a 
noisy (possibly very bad) estimate of the channel is available at the 
receiver. First, a decoding metric based on the channel a posteriori 
probability density, conditioned on the channel estimate is derived 
and used for decoding BICM multiband OFDM. Then, we charac- 
terize the limits of reliable information rates in terms of the maxi- 
mal achievable outage rates associated to the proposed metric. We 
also compare our results with the outage rates of a system using a 
theoretical decoder. Our results are useful for designing a communi- 
cation system where a prescribed quality of service (QoS), in terms 
of achievable target rates with small error probability, must be satis- 
fied even in the presence of imperfect channel estimation. Numeri- 
cal results over both realistic UWB and theoretical Rayleigh fading 
channels show that the proposed method provides significant gain in 
terms of BER and outage rates compared to the classical mismatched 
detector, without introducing any additional complexity. 

1. INTRODUCTION 

Ultra- Wide-Band (UWB) is defined as any wireless transmission 
scheme that occupies a bandwidth of more than 25 % of its cen- 
ter frequency or greater than 500 MHz over the 3.1-10.6 GHz fre- 
quency band | Tj. Multiband Orthogonal Frequency division multi- 
plexing (MB-OFDM) [2 1 is a spectrally efficient technique proposed 
for high data rate, short range UWB applications. This approach 
uses a conventional OFDM system, combined with bit interleaved 
coded modulation (BICM) and frequency hopping for improved di- 
versity and multiple access. In MB-OFDM, the channel is assumed 
changing so slowly that it is considered time invariant during the 
transmission of an entire frame. Channel estimation is performed by 
one known symbol (pilot) transmitted at the beginning of the infor- 
mation frame while the rest of the frame is decoded based on the 
estimated channel. Due to the limited number of pilots, the estimate 
of the channel is imperfect and the receiver has only access to this 
noisy channel estimate. However, the receiver/decoder metric for 
any maximum-likelihood (ML) based detector, requires knowledge 
of the exact channel. 

A standard sub-optimal technique, known as mismatched ML 
decoding, consists in replacing the exact channel by its estimate 
in the receiver metric. Hence, the resulting decoding metric is not 
adapted to the presence of channel estimation errors (CEE). This 
practical scenario, with mismatched decoding, leads to the following 
important questions. Firstly, in presence of CEE with a given amount 



of training, what are the limits of reliable transmission (the capacity). 
Secondly, what type of encoder/decoder is necessary to transmit re- 
liable information close to the performance limits. The first problem 
has been recently addressed in [3], characterizing the maximal reli- 
able information rate, by using the notion of estimation induced out- 
age capacity. Unfortunately, the theoretical encoder/decoder used to 
achieve this capacity can not be implemented on practical commu- 
nication systems. Besides, the mismatched ML decoding has been 
showed to be largely suboptimal for the considered class of chan- 
nels. Basically, the transmitter and the receiver strive to construct 
codes for ensuring reliable communication with a quality of service 
(QoS), no matter what degree of channel estimation accuracy arises 
during the transmission. The QoS requirements stand for achieving 
target rates with small error probability even with very bad chan- 
nel estimates. In this paper, we propose a practical decoding met- 
ric for the aforementioned mismatch scenario. This metric uses the 
estimated channel and the a posteriori pdf characterizing the chan- 
nel estimation process, which matches well the channel knowledge 
available at the receiver. Based on the derived metric, we formulate 
our decoding rule for BICM MB-OFDM. Interestingly, the present 
metric coincides with that derived for space-time decoding from in- 
dependent results in |4). In order to determine the limits of reli- 
able information rates associated to the proposed metric, we use the 
complementary results obtained in |5|. This allows us to compare 
the maximal supported rate associated to our metric with that of the 
classical mismatched ML and the theoretical decoder. Our results 
are relevant for communication systems where a prescribed quality 
of service (given by the outage probability) must be ensured even in 
the presence of imperfect channel estimation. 

The outline of this paper is as follows. In Section|2]we describe 
the system model for MB-OFDM transmission over a frequency se- 
lective fading channel. Section[3]presents the pilot assisted channel 
estimation: we specify the statistics of the CEE and then calculate 
the posterior distribution of the perfect channel conditioned on the 
estimated channel. This posterior distribution is used in section|4]to 
derive the ML decoding metric in the presence of imperfect channel 
state information at the receiver (CSIR). In section|5] we use the gen- 
eral modified metric for soft decoding BICM MB-OFDM systems. 
In section [6] we derive the achievable outage rates of a receiver us- 
ing the proposed metric. Section |7] illustrates via simulations the 
performance of the proposed receiver in realistic UWB channel en- 
vironments and section[8]concludes the paper. 

Notational conventions are as follows : upper case bold symbols 
denote vectors, Ijv represents an (iV x iV) identity matrix; (.) T and 
denote vector transpose and Hermitian transpose, respectively. 
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Fig. 1. TX architecture of the multiband OFDM system. 

2. TRANSMISSION MODEL 

An MB-OFDM system divides the spectrum between 3.1 to 10.6 
GHz into several non-overlapping subbands each one occupying ap- 
proximately 500 MHZ of bandwidth |2|. Information is transmitted 
using OFDM modulation over one of the subbands in a particular 
time-slot. The transmitter architecture is depicted in figure [T] each 
bit bi is convolutionally encoded into two bits c° and c\ which are 
interleaved in order to break the error bursts. The interleaved bits 
are gathered in subsequences of B bits dj, ■ ■ ■ , d*. and mapped to 
complex M-QAIVQ (M = 2 B ) symbols S k with average energy 
E s — E[\S k \ 2 ]. At a particular time-slot, a time-frequency code 
(TFC) selects the center frequency of the subband over which the 
OFDM symbol is sent. The TFC is used not only to provide fre- 
quency diversity but also to distinguish between multiple users. 

At the receiver, assuming a cyclic prefix (CP) longer than the 
channel maximum delay spread and after performing fast Fourier 
transform (FFT), OFDM converts the the channel into M parallel 
Rayleigh distributed flat fading subchannels, where M is the number 
of subcarriers. The baseband equivalent observation model over a 
given subband can be written as 

Y = D H S + Z, (1) 

where (M x 1) vectors Y and S denote received and transmitted 
symbols, respectively; the noise block Z is assumed to be a circularly 
symmetric complex Gaussian random vector with distribution Z ~ 
CM(0, a 2 1m ) ; and D H is a diagonal matrix with diagonal elements 
equal to H = [Hi, . . . , Hn] T , where H is the vector of channel 
FFT coefficients. 

3. PILOT ASSISTED CHANNEL ESTIMATION 

Under the assumption of time-invariant channel characteristics over 
the entire transmitted frame, channel estimation is usually performed 
by transmitting known training symbols at the beginning of each 
frame. 

Assume that we are interested to estimate the fc-th, k = {1, 
M}, fade coefficient Hk via the transmission of N pilot symbols 
and let S T = [Si, ... , Sn] T be the (N x 1) vector of transmitted 
training sequence. We will assume constant modulus training sym- 
bols for the fc-th subcarrier, i.e., \Si | 2 = ... = |<S]v| 2 = Pr,k- From 
we have 

Y T =D Hfc S T +Z T , (2) 
where D Hfc is a diagonal matrix with diagonal elements equal to the 
(iV x 1) vector H k = [Hk, ■ ■ ■ , Hk] T and Z T is a noise vector 

'Here, 16-QAM mapping is used instead of QPSK, proposed in |6|. 



with the same distribution as Z affecting the transmission of pilot 
symbols. The ML estimate of Hk is obtained by maximizing the 
likelihood function p (Y T \Hk, S T ). After some standard calculus 
we have 

m iL = ^I^=Hk + £k, fc = l,...,M (3) 

JV r T ,k 

where 8k ~ CM(0, a% ) is the estimation error. From equations l[3), 
it is obvious that cr| fc = SNR"^, where SNR T , fc = 

Assuming that HkS are i.i.d. and distributed as CM(0, a 2 ), the 
posterior distribution of the perfect channel conditioned on the esti- 
mated channel coefficient is given by p(Hk\Hk) oc p(Hk\Hk)p(Hk) 
After standard manipulation of Gaussian densities we obtain 

p(Hk\H k ) = CM (^pHk , . (4) 

where p k = o\l{ai + a% k ). 

The availability of the estimation error distribution constitutes 
an interesting feature of pilot assisted channel estimation that we 
used to derive the posterior distribution ©. Next, we see how this 
additional information is exploited in a modified ML metric for im- 
proving the detection process. 

4. MAXIMUM-LIKELIHOOD DETECTION METRIC IN 
THE PRESENCE OF CHANNEL ESTIMATION ERRORS 

4.1. Mismatched ML metric 

The observation model Q} can be rewritten in a component-wise 
form as 

Y h = H h S k + Z h , k = l,...,M (5) 

It is well known that the ML estimate of Sk under the assumption 
of Gaussian noise is given by maximizing p(Yk\Hk, Sk) which is 
equivalent to minimizing the euclidean distance 

Sf L (Hk) = arg min { \Y k -H k S k \ 2 }, (6) 

where the set 5? contains all of the possible discrete values in the 
constellation that Sk can take. Note that the above metric depends 
on the realization of the perfect channel Hk and is optimum when 
perfect CSIR is available. However, in a real communication system, 
the receiver has only access to an imperfect estimated version of the 
channel. In the mismatched decoder, the estimated channel Hk is 
used instead of Hk in the ML metric of ([SJ as 

Sf L (i? fc ) = argmin {\Y k - H k S k \ 2 } \ , (7) 

which leads to a sub-optimal solution due to the mismatch intro- 
duced by the CEE. Next, we will see how the derived posterior chan- 
nel distribution I0, can be used in a modified ML metric that takes 
into account the available imperfect CSI. 

4.2. Modified ML metric for imperfect CSIR 

Since in practice the receiver knows solely the imperfect channel, 
it is relevant to express the likelihood criterion in terms of the esti- 
mated fade coefficient. For this purpose, we consider the problem of 



detecting the transmitted symbol S k by using the imperfect channel 
estimate H k = H ky o in an ML criterion as 

S^(H k , ) = argmin V M (Y k ,S k \H k = H k , ), (8) 

V/hereV M (Y k ,S k \H k ) 4 -kip(Y k \H k , S k ). Thep<tfp(Y k \H k , S k ) 
can be calculated as 



p(Y k \H k ,S k ) = f p(Y k ,H k \H k ,S k ) dH k , 

J Hi. 



(9) 



The joint likelihood function p(Y k , H k \H k ,S k ) is related to known 
probabilities as 

p(Y k ,H k \H k ,S k )=p{Y k \H k ,S k )p(H k \H k ), (10) 

where the last equation results from the independence between S k 
and (H k ,H k ). It is clear from ® and ([Toj that p(Y k \H k ,S k ) is 
obtained by the following conditional expectation 

p(Y k \H k ,S k )=E HklBk [p(Y k \H k ,S k )\H k ], (11) 

evaluated over the posterior distribution of the exact channel given 
its estimate (equation $4$). Since p(H k \H k ) and p(Y k \H k , S k ) are 
Gaussian densities, their product remains Gaussian and it is easy to 
verify that p(Y k \H k , S k ) = CN{p M ,(J 2 M ) where p M = pH k S k 
mdo 2 M =o 2 + (l-p)\S k \ 2 . 

According to <[8j, the decision metric to be minimized for ML 
detection at the receiver is easily seen to be 



V M {S k ,Y k \H k ) =\n{<jl + {l-p)\S k \ 2 



\Y k 



pH k S k \ 2 



a'i + (l- p)\S k \ 2 ' 
(12) 

Let us first consider the ideal situation when exact channel knowl- 
edge is available. In this case p{H k \H k ) reduces to S(H k — H k ) 
and the averaging E H| _ (.) of i ll It becomes equivalent to replacing 

H k by H k in p(Y k \H k , S k ). Similarly, we note that with near per- 
fect channel (i.e. <r| — > or TV — > do), p — > 1 and we have 
T> M (Y k , S k \H k ) — > \Y k — H k S k \ 2 , which is the euclidean dis- 
tance metric for perfect CSIR. However, in the presence of CEE, 
the modified metric evaluates the average of the pdf p(Y k \H k , S k ) 
(that would be used if the exact channel were known) over the true 
channel H k using the conditional expectation in i ll It . 

While we have focused on OFDM, the metric l !12t can also be 
applied for ML decoding/detection of single carrier systems with fre- 
quency selective fading. In the following, we use the above modified 
metric in the decoding process of BICM MB-OFDM. 

5. DECODING OF BICM-OFDM BASED ON IMPERFECT 
CHANNEL ESTIMATION 

The problem of decoding BICM-OFDM with large size constella- 
tions has been addressed in 1 7 1 under the assumption of perfect CSIR. 
Here we consider the same problem in the case where the decoder 
has an imperfect estimate of the channel state. We denote by S\ = 
[Sf , . . . , S^] the sequence of transmitted complex symbols and by 
Yi = [Yf , . . . , Y^] the corresponding sequence of received sym- 
bols where r is number of OFDM symbols in a frame. Based on the 
transmission model of figure Q] the sequence of transmitted bits b 
has to be estimated in the ML sense so as to maximize the likelihood 



function p(Yl b). In practice, due to the one to one correspon- 
dence between the information bits and the interleaved codeword d, 
one look for the ML estimate d of the interleaved codeword as 



argmaxp(YI|HI,SI). 



(13) 



By assuming that the received symbols are independent and the chan- 
nel noise is white, U3\ is equivalent to (see [7|) 

B 

d M « argmaxT] T\p(d l k \H k ,Y k ), (14) 
where the bit metrics p(d l k \H k , Y k ) are given by 



p(d{\H k ,Y k ) oc P(Y k \H k ,S k (dl,...,d%)), 
d£e{o,i} 



exp I 



(15) 



dJe{o,i} M 



and p(Y k \H k , S k ) is specified from dill ). We emphasize that the 
above decoding rule is conditioned on the imperfect estimated chan- 
nel sequence H[ available at the receiver. Information bits are es- 
timated in the trellis decoder part from the above bit metrics by ap- 
plying either non-iterative or turbo procedures (see figure^. 
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Fig. 2. BICM iterative demodulation with CEE. 



6. ACHIEVABLE OUTAGE RATES ASSOCIATED TO THE 
PROPOSED METRIC 

We now derive the achievable rates C M of a receiver using the 
modified metric T> M of j 1 2b . This is done by using the follow- 
ing theorem proposed in [5 ] for a discrete-time memoryless channel 
W(y|s,H) = nf =1 CAA(ff fcSfc ,a z 2 ), wherey = (y u ...,y M ), 
s = (si, . . . , sm) with decoding metrics T> M (y k , s k \H k ): 

Theorem 1. 

C 



, M - sup inf D(foP s \\P s P Y ) 



where the relative entropy functional D(JoPs\\PsPy) = 

/(y|s) 



// 



P s (s)/(y|s)log- 



-dyds, 



(16) 



(17) 



JPs'(8')f(y\B')daf 
and T denotes the set of all possible channels /(y |s) that satisfy: 
a : E s [/(y|s)l =K.[W(y|B,H)l, a.s. (18) 



<-j : JK /(y|s) [D(s k ,y k \H k ) 
for every k — 1 , . . . , M, 



E w (y| s) [V(sk,y k \H k ) 
(19) 



We assume Gaussian i.i.d. inputs Sk ~ CJ\f(0,P). Further- 
more, since we aim to minimize the mutual information in l |16t , we 
choose the Gaussian distribution /(y|s, (j) = Ylk=i CM{pkSk, 0fc) 
with n — (fix, . . . , (j,m)- It is easy to see from constraint ci that the 
variances of the optimal channel distribution /(y|s, p) are given by 



16-QAM, 3 decoding iterations, channel CM1 



= o" z + P(\H k \ -\nk\), k = l, 



,M. 



(20) 



Then, for any given channel H and its estimate H, the resulting 
achievable rates for the proposed metric are obtained as 



C M {Ti,&)= inf _ £log 2 (l + 



(21) 



where V(Hk, Hk) are the constraint sets C2 corresponding to the 
metric T>m , derived from the inequality constraint l |19t as 

V(H k ,Hk) = {fik e C:\Hk- a k H k \ 2 < \ft k - aH k \ 2 }, (22) 

where a k = pk{^ko'i~ P(l — Pk) )/(Afc<7 2 -P(l — Pfc)(l — Afc) ); 
A, = exp (^2- T )£ 1 (^L- T ) and ^(s) 4 /+« -ipOdu 
is defined as the exponential integral. 

We note that j_ 1 | i is an increasing function of ||/x|| 2 . Conse- 
quently, the infimum is obtained for fik = p, op t,k resulting from 
minimizing |/i|| 2 constrained by the equality in the constraint set 
l !22t . Finally, by using Lagrange multipliers we obtain 



k = 1, 



where 



= < 



a k + 



\H k — a.kHk\ 

\H k \ ^ 
\H k — a.kHk\ 

\Hk\ 
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if ak > 
if a k < 0. 



(23) 



(24) 



The achievable rates associated to the metric T> M are given by 

M / p 2 I rr |2 

C M (H,H) = Vlog 2 1 + fcl ^ 

(25) 

Then, the associated outage probability for an outage rate R > is 
defined as 

P^ ut (i?,H) = Pr H| H(H G A^HJIH), 

with (i?, H) = {H : C M (H, H) < ii} and the maximal 
outage rate for an outage probability 7 is 

C^ ut (7,H) = sup{i? > : P^(R,&) < 7}. (26) 

Since this outage rate still depends on the random channel estimates 
H, we will consider the average rates over all channel estimates 

C(7)=%[C(7-H)]. (27) 

The achievable rates i27\ are upper bounded by the estimation-induced 
outage capacity, the maximal outage rates (i.e. maximizing over all 
possible transmitter-receiver pairs), derived in [3] by using a theo- 
retical decoder. In our case, this capacity is given by 



Cr(7)=E fi [C° ut (7,H)], 
where the outage rates Cg U (7, H) are computed using 



C s (H,H)=]Tlog 2 1 + 



P\Hk\ 2 



(28) 



(29) 
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Fig. 3. BER performance of the proposed decoder over the LOS 
UWB channel CM1 for various training sequence lengths. 
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Fig. 4. BER performance of the proposed decoder over Rayleigh 
fading channel for various training sequence lengths. 



7. SIMULATION RESULTS 

Simulations have been carried out in the context of IEEE802.15 
wireless PAN [6|: M — 100 data subcarriers along with 32 CP 
samples compose one OFDM symbol. Information bits are encoded 
by a rate R = 1/2 convolutional encoder with constraint length 3 
defined in octal form by (5,7). The interleaver is a pseudo-random 
one operating over the entire frame and the mapping is 16-QAM. 
For each frame that contains 100 OFDM symbols, a different real- 
ization of the realistic UWB channel model specified in 1 8 1 has been 
drawn and supposed constant during the whole frame. Besides, the 
same average energy is allocated to training symbols and data, i.e., 
E T = E(\S k \ 2 ) = E s . 

Figure [3] depicts the bit error rate (BER) performance gain that 
is obtained by decoding BICM MB-OFDM with the modified ML 
decoder in the context of LOS CM1 channel, where the channel es- 
timation is performed by sending different training sequences with 
lengths N = {1, 2, 8} per frame. It can be noticed that the proposed 



decoder outperforms the mismatched decoder especially when few 
numbers of pilot symbols are dedicated for channel estimation. Note 
that the modified decoder with TV — 1 pilot performs very close to 
the mismatched decoder with TV = 2 pilots. For comparison, results 
obtained with theoretical Rayleigh fading channel are illustrated in 
figure|4] It can be observed that for TV = 2, the SNR to obtain a BER 
of 10 -3 is reduced by about 1.5 dB if the modified ML decoding is 
used instead of the mismatched approach. 

Figure [5] shows the BER performance versus the training se- 
quence length TV at a fixed Eb/N of 12 dB for the CM1 channel. 
This allows to evaluate the number of training sequence necessary 
to achieve a certain BER. We observe that at Eb/No = 12 dB, the 
modified ML decoder requires TV = 9 pilot symbols per frame to 
achieve a BER of 1CP 4 while the mismatched decoder attains this 
BER for TV = 12 pilot symbols. Besides, this has been outlined 
in Section IV, for large training sequence lengths (TV > 12), both 
decoders have close performance. 

Figure[6]shows the expected achievable rates j27\ corresponding 
to one OFDM symbol transmission with M = 16 data subcarriers 
(in bits per channel use) versus the SNR, for an outage probability 
7 = 0.01. For comparison, we show the expected outage rates of 
mismatched ML decoding and that obtained in the case of the perfect 
theoretical encoder/decoder J28l >. It can be observed that at a mean 
outage rate of 4 bits, the modified metric reduces the SNR by about 1 
dB. Indeed, we observe that the achievable rates using the proposed 
decoding metric are close (only 2 dB of SNR) to those obtained by 
using a theoretical decoder. 

8. CONCLUSION 

This paper studied the problem of ML reception when the receiver 
has only access to a noisy estimate of the channel in the case of pilot 
assisted channel estimation. By using the statistics of the estima- 
tion error, we proposed a modified ML criterion that is expressed in 
terms of the estimated channel coefficient. This modified metric let 
us to formulate appropriate branch metrics for decoding BICM MB- 
OFDM with imperfect channel knowledge. Achievable outage rates 
associated to this modified metric have also been derived. Compar- 
ison in terms of BER and achievable outage rates was made with 
a system whose receiver uses mismatched ML decoding and a per- 
fect theoretical decoder. Simulation results conducted over realis- 
tic UWB channels, indicate that mismatched decoding is quite sub- 
optimal for short training sequence and confirmed the adequacy of 
the proposed decoding rule in the presence of channel estimation 
errors. This was obtained without introducing any additional com- 
plexity. 
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Fig. 5. BER versus the number of pilot symbols at Eb/No = 12 dB 
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